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Abstract 

We used an ITS2 primary and secondary structure and Compensatory Base Clianges (CBCs) analyses on new 
French and Spanish Dunollela solino strains to investigate their phylogenetic position and taxonomic status within 
the genus Dunoliello. Our analyses show a great diversity within D. solino (with only some clades not statistically 
supported) and reveal considerable genetic diversity and structure within Dunoliello, although the CBC analysis did 
not bolster the existence of different biological groups within this taxon. The ITS2 sequences of the new Spanish 
and French D. solino strains were very similar except for two of them: ITC5105 "Janubio" from Spain and ITC51 19 
from France. Although the Spanish one had a unique ITS2 sequence profile and the phylogenetic tree indicates 
that this strain can represent a new species, this hypothesis was not confirmed by CBCs, and clarification of its 
taxonomic status requires further investigation with new data. Overall, the use of CBCs to define species 
boundaries within Dunoliello was not conclusive in some cases, and the ITS2 region does not contain a 
geographical signal overall. 

Keywords: Canary Islands, Compensatory Base Changes, Dunoliello solino, Internal Transcribed Spacer, Saltworks, 
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Background 

The Internal Transcribed Spacer 2 (ITS2) of the nuclear 
rDNA cistron is one of the most frequently used regions 
for phylogenetic analysis in algae [1-3]. Although its 
application in deep taxonomic levels was initially limited 
to comparisons of genera within the same family owing 
to uncertainties in alignment at higher taxonomic levels, 
the analysis of its secondary structure has provided key 
solutions to this problem [4]. Thus, the use of an ITS2 
secondary structure improves sequence alignments, 
resulting in a higher robustness and accuracy of phylo- 
genetic reconstructions [5] and providing help to distin- 
guish species [6]. Furthermore, an automatic approach 
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to analysis is possible [7], as a pipeline consisting of the 
ITS2 Database (annotation/structure prediction), 4SALE 
(alignment), ProfDistS (inferring phylogenies) and the 
CBCAnalyzer (distinguishing species) have recently 
become available (http://its2,bioapps.biozentrum,uni- 
wuerzburg.de/?about). 

In Dunaliella (Chlorophyceae), the use of ITS2 sec- 
ondary structure for phylogenetic analysis has a long 
tradition [8-13]. The genus Dunaliella comprises 
twenty-eight recognized species separated in two subge- 
nera, Pascheria (which contains the freshwater species), 
and Dunaliella (grouping the marine species); the latter 
is further subdivided into four sections: Tertiolecta, 
Dunaliella, Viridis and Peirceinae [11]. The species 
ascribed to these four sections occur in a wide range of 
marine habitats such as oceans, brine lakes, salt 
marshes, salt lagoons and salt water dishes near the sea 
[14], being Dunaliella salina Teodoresco (section 
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Dunaliella) the most studied one. Dunaliella salina, is 
the most halotolerant eukaryotic photosynthetic organ- 
ism known to date [14,15] since it shows a remarkable 
degree of adaptation to a variety of salt concentrations 
and it accumulates large amounts of carotenes under 
extremely stressful conditions such as high salinity, low 
nitrogen levels, and high solar radiation [14]. Nowadays, 
it is the best commercial source of natural p -carotene 
[14,16], and it also stands out as a source of glycerol 
[17]. 

One of the aspects of D, salina that have most intri- 
gued researchers is the enormous variability within 
strains regarding its geographic, physiological, and mor- 
phological characteristics [18-24]. Recent phylogenetic 
analyses of ITS1+ITS2 combined with the analysis of 
the ITS2 secondary structure of D, salina strains have 
also revealed a high intraspecific variation [8-11,25]. 

The high genetic diversity detected in D, salina, its 
morphological plasticity, and the restricted geographical 
sampling used in all scientific publications to date have 
hindered taxonomic elucidation in this taxon. Our 
objectives are (1) to use the ITS2 sequences and second- 
ary structure analysis in a thorough geographic and 
taxonomic representation of the genus Dunaliella and 
particularly D, salina (including new Spanish and 
French strains) to improve understanding of the com- 
plex phylogenetic structure in this taxon, (2) to study 
the relationship of the new strains with the Dunaliella 
sequences available at the ITS2 Database and/or at Gen- 
Bank [26]; (3) to investigate if the Compensatory Base 
Changes (CBCs) analyses could elucidate the species 
concept in Dunaliella, and reveal potentially new spe- 
cies; and (4) to test if D, salina ITC5105 "Janubio" 
could be considered a new species. 

Results 

The ITS2 primary and secondary structure phylogenetic 
analysis of all the Dunaliella sequences available in the 
ITS2 Database plus the new sequences revealed great 
heterogeneity, although some of the clades were not sta- 
tistically supported (Figure 1). No phylogenetic relation- 
ship is supported between the two Dunaliella freshwater 
species, since D. lateralis was positioned outside the 
Dunaliella subgenus, while D, acidophila was positioned 
within this subgenus. 

The strains positioned in the tertiolecta-clade repre- 
sent species that belong to different traditionally 
accepted sections: section Tertiolecta [D, tertiolecta, D, 
primolecta, D, quartolecta, D. polymorpha), section Viri- 
dis {D. bioculata, D, minuta), section Dunaliella [D, sal- 
ina, D. parva), and section Peircei {D. percei). The 
majority of these strains had an exact ITS2 sequence 
(Figure 1). Most D. viridis strains sampled were posi- 
tioned in a single clade; however, some strains 



(UTEX200, CCAP19/3) were positioned in a different 
clade together with D, pseudosalina CONCOlO (pseudo- 
salina-clade), but without statistical support. 

The majority of the D. salina strains were distributed 
in two different clades (salina-clade-I and salina-clade 
II), positioned between the tertiolecta sub-clades, 
whereas two D. salina strains (CCAP19/30, CCAP19/18) 
were positioned together in a different clade. Only two 
of the new ITS2 D. salina sequences, ITC5119 ("salina- 
clade-I") and D. salina ITC5105 "Janubio" (salina-clade 
II), had a unique ITS2 sequence profile (Figure 1). The 
CBC analysis of the ITS2 secondary structure showed 
that there is at least one CBC between "Janubio" and 
the other D, salina strains analyzed (see Additional file 
1), except D. salina CCAP19/30 and CCAP19/18; how- 
ever these two strains were phylogenetically distant to 
Janubio. 

The taxonomic identification of the Dunaliella sp. 
strains (below the D. acidophila strain in Figure 1) was 
not possible because they did not match with any 
known Dunaliella species. Despite the phylogenetic tree 
indicates that this group may indeed represent a new 
species; the analysis of the CBC did not confirm these 
results (Additional file 1). 

We observed a lack of CBCs in some strains when we 
compared species in different clades: viridis vs salina, 
viridis vs tertiolecta, salinas vs tertiolecta, salina vs sal- 
ina, etc (Additional File 1). Also, the CBC analysis of the 
different species within the subgenus Dunaliella for the 
confirmation of the species boundaries was not conclu- 
sive in some cases (Additional file 1). 

The phylogenetic results of the ITS2 sequences in this 
study, and previous confirmed identification of some 
Dunaliella strains (see Table 1), allow us to suggest the 
re-identification of some strains (Table 1). 

Finally, the ITS2 data was not informative regarding 
the geographic origin of the D. salina strains. 

Discusion 

The ITS2 Database allows automatic large scale simulta- 
neous analyses of both ITS2 sequences and their sec- 
ondary structures. Potential pitfalls are in structures 
obtained by different algorithms; however the main diffi- 
culty of performing a phylogenetic analysis of the genus 
Dunaliella is the misinformation available at Culture 
Collections and GenBank regarding the identification of 
strains and sequences. The tracking of the true identifi- 
cation of each strain was only possible after consulting 
several publications where the authors concluded that 
they were misidentified and suggested their re-identifica- 
tion [9,10,27]. In this study, we have also suggested the 
re-identification of some strains based on our ITS2 data. 
To overcome all these unnecessary problems (and given 
that it is unlikely that all the Dunaliella strains could be 
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I Oedogoniales (outgroup) 



71482662 Oedogonium nodulosum 

- 62183531 Bulbochaete rectangularis var. Hiloensis UTEX954 

- 2982755 Astrephomene gubernaculifera UTEX2479 
■ 14165185 Pandorina morum Ep3 




111073405 Chlamydomonas reinhardtii SAG11-31 
- 4007503 EudorinaelegansUTEX1192 
■ 14091689 Yamagishiella unicocca UTEX2431 
18025134 Phacotus lenticularis 970805-20 
2627275 Lobochlamys culleus UTEX1057 
16596847 Dunaliella lateralis Nepal •— 



J 87047580 Dunalif 



"salina-clade I" 



tertiolecta-clade 



"salina-clade M' 



Pseudosalina-clade 



13274452 Chlamydomonas allensworthii Isolate 7 
87047578 Dunaliella sp. CCMP367 

197290646 Dunaliella sp. ABRIINW Ml/2 

Dunaliella acidophila CCAP19/35 
87047594 Dunaliella sp. SPMO 112-3* 

87047601 Dunaliella sp. SPMO 201-3 i — = 

87047596 Dunaliella sp. SPMO 128-2 • rtS^ n""1'1!^ 

"100 1 87047591 Dunaliella sp. SPMO 109-1 87047603 Dunaiteiia Id 

' 87047586 Dunaliella sp. CCMP220 I 

87047604 Dunaliella sp. SPMO 201-6 
87047590 Dunaliella sp. FL1 

87047592 Dunaliella sp. SPMO 112-1 87047593 Dunaliella s 

87047609 Dunaliella sp. SPMO 300-4 
87047607 Dunaliella sp. SPMO 210-3 



.a 



IV 



SPMO 200-3 
SPMO 201-4 
SPMO 201-5 



^ 1 



205361369 Dunaliella salina CCAP19/30 •— 
145587829 Dunaliella salina CCAP19/18 
87047599 Dunaliella sp. SPMO 200-8 
87047589 Dunaliella viridis CONC002 
205371719 Dunaliella sp. ABRIINW U2/1 
J- 87047612 Dunaliella sp. SPMO 601-1 

P— 87047597 Dunaliella sp. SPMO 200-2 • 

I — J- 61200914 Dunaliella viridis SAG44.89 
Uif— 58339343 Dunaliella viridis SHU 
~i| — 56578597 Dunaliella viridis DCCBC4 
^^"^i- 61200913 Dunaliella viridis DCCBC3 

4 16596845 Dunaliella parva UTEX1 983 

I 87047572 Dunaliella sp. BSF1 • 

n_r 87047605 Dunaliella sp. SPMO 202-4 •— 

^ 16596834 Dunaliella viridis CONC002 

Dunaliella salina ITC5119 

2645739 Dunaliella sp. 006 
16596837 Dunaliella salina CONC006 
16596840 Dunaliella bardawil ATCC30861 
47499297 Dunaliella salina DCCBC1 
Dunaliella salina CCAP19/12 « 



-\ 33333776 Dunaliella s£ 



U-U mismatch I ^ 

ii *«». 



87047574 Dunaliella s| 

87047610 Dunaliella s| 

87047611 Dunaliella s( 
. 87047573 Dunaliella s| 



BSF3 

SPMO 300-5 
SPMO 600-1 
BSF2 




71482594 Dunaliella salina OUC36 
33333772 Dunaliella salina Yucatan 
33333774 Dunaliella salina Tanggu •— 

Dunaliella percei [tertiolecta] CCAP19/2 
213958821 Dunaliella salina DS18S1 

145587823 Dunaliella tertiolecta Dtsi • 

Dunaliella bioculata CCAP19/4 •— 



145587824 Dunaliella salina Dsge 

145587825 Dunaliella tertiolecta SAG13.86 

145587826 Dunaliella salina SAG42.88 

145587827 Dunaliella tertiolecta ATCC30929 
145587830 Dunaliella tertiolecta CCAP19/27 
87047577 Dunaliella bioculata UTEX199 

87047583 Dunaliella primolecta UTEX1000 

87047584 Dunaliella sp. CCMP1641 

87047585 Dunaliella sp. CCMP1923 

87047587 Dunaliella tertiolecta CCMP1302 

87047588 Dunaliella tertiolecta CCMP364 



minutaCCAP19/S 
tertiolecta CCAP1 9/23 
Dunaliella tertiolecta CCAP19/6B 



16596844 Dunaliella parva [tertiolecta] CCAP19/9 

16596842 Dunaliella tertiolecta UTEX999 

16596841 Dunaliella tertiolecta CCMP1 320 • 

2627284 Dunaliella tertiolecta FHL 
56578596 Dunaliella tertiolecta DCCBC5 
213958823 Dunaliella salina DS18S3 
213958822 Dunaliella salina DS18S2 

Dunaliella salina ITC5105 "Janubio" 
47499296 Dunaliella salina AC144 

■ 51 573027 Dunaliella salina 9802 Dunaliella salina itc5118 

16596839 Dunaliella salina UTEX1644 • Dunaliella salina itc5122 

33333770 Dunaliella salina CONC001 Dunaliella salina ccapi9/25 

16596838 Dunaliella salina CONC007 
16596835 Dunaliella pseudosalina CONC010 
197290927 Dunaliella sp. ABRIINW Ml/1 
|_f- 145587828 Dunaliella salina [viridis] CCAP19/3 
\r- 87047582 Dunaliella parva [maritima] SAG19-1 
li- 49073091 Dunaliella salina SAG184.80 
JJ' 16596836 Dunaliella salina [viridis] UTEX200 



254838316 Dunaliella salina MSI-1 

71482595 Dunaliella salina OUC38 

71482596 Dunaliella salina OUC66 

71482597 Dunaliella salina hd5 

71482598 Dunaliella salina hd6 
205371718 Dunaliella sp. ABRIINW U/1 
46250926 Dunaliella salina DCCBC2 
71482593 Dunaliella salina 0UC21 
Dunaliella salina ITC51 03 
Dunaliella salina ITC5000 
Dunaliella salina ITC5003 
Dunaliella salina ITC5100 
Dunaliella salina ITCS101 
Dunaliella salina ITCS102 
Dunaliella salina ITC5104 
Dunaliella salina ITC5106 
Dunaliella salina ITC5107 
87047576 Dunaliella bardawil UTEX2538 



^5l 



16596843 Dunaliella parva [tertiolecta] CCMP362 
16596846 Dunaliella pelrcel [tertiolecta] UTEX2192 
47933783 Dunaliella tertiolecta CCAP19/6B 
50952902 Dunaliella primolecta UTEX1000 

51035302 Dunaliella minuta SAG23.86 

51035303 Dunaliella maritima SAG42.89 

51035304 Dunaliella sp. SAG19.6 
55979209 Dunaliella tertiolecta CCAP19/27 

71482599 Dunaliella bardawil hd7 

71482600 Dunaliella primolecta hd8 

71482601 Dunaliella parva hd9 

71482602 Dunaliella sp. hdIO 
76097092 Dunaliella bioculata UTEX199 

77539931 Dunaliella polymorpha CCAP19/7C 

77539932 Dunaliella quartolecta CCAP19/8 
77955899 Dunaliella tertiolecta DCCBC26 



Figure 1 Sequence-structure tree and consensus structure of ITS2. (A) Sequence-structure Neighbor-Joining tree obtained by ProfDistS and 

supporting bootstrap values (100 replicates). Strains on the right written within squares have identical ITS2 sequences. Strains written in bold 

were sequenced in this study. All strains information can be found in Supplementary Table 1 (B) Consensus structure (51%) for all ITS2 

sequences obtained from the complete multiple sequence-structure alignment without gaps. Helices are numbered l-IV. Sequence conservation 

is indicated from red (not conserved) to green (conserved). Nucleotides which are 100% conserved in all sequences are written as A, U, G or C. 

Nucleotide bonds which are 100% conserved throughout the alignment are marked in yellow. The figure was generated with 4SALE [41]. 
^ J 



openly available), we suggest to establish a "type strain" 
for each Dunaliella taxon (including subspecies, forms 
or varieties). These basic data should be easily obtained 
from any official culture collection, thereby greatly facili- 
tating comparison with new field isolates and avoiding 
misleading information and/or false conclusions. 

Our ITS2 phylogenetic analysis of Dunaliella reveals 
several major groups, and positions the freshwater D, 
lateralis clearly outside Dunaliella, confirming that it no 
longer should be considered a member of this genus 
[10]. Nevertheless, the other freshwater species analyzed 
in this study [D, acidophila, CCAP19/35), maintained its 
position within the subgenus Dunaliella, and was not 
phylogenetically related to D, lateralis, as recently pro- 
posed [28]. Furthermore, the observation that different 
species belonging to several Sections {Tertiolecta, Viri- 
dis, Dunaliella and Peircei) share the exact same ITS2 
sequence, make us believe that they correspond to a sin- 
gle species. These data agree with other authors 
[10,11,27,29], who suggested that the number of 



Dunaliella species may be much lower than it has been 
claimed till now. The possibiUty that the ITS2 gene is 
not able to discriminate between these species is highly 
unlikely; therefore, our observations support the sugges- 
tion that the morphological and physiological criteria 
available to discriminate Dunaliella species are either 
not very reliable [11], or are difficult to interpret. 

In an attempt to clarify the species concept within 
Dunaliella, we searched for compensatory base changes 
(CBCs). Several case studies have revealed that the 
detection of a CBC in the ITS2 secondary structure 
between two organisms is correlated with sexual incom- 
patibility [6,30,31], and these changes have been pro- 
posed as markers for distinguishing species [6,7,30,31]. 
In summary, these investigations conclude that while a 
CBC in a pair of sequences is positively correlated with 
species distinctness at a confidence level of 93%, the 
lack of a CBC in the ITS2 secondary structure does not 
necessarily indicate that two organisms belong to the 
same species [6]. The overall analysis of the CBC was 
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Table 1 Taxonomic classification. Culture Collection, Geographic Origin and GenBank acession numbers of the strains 
included in this study. 



Former classification 



New suggested classification 
[Reference] and comments 



Culture 
Collection 



Geographic Origin 



Isolator GenBank GenBank 
(Date) acession identifier 

number 



Subgenus Poscherio (Freshwater species) 

Dunoliello lateralis 
Pascher SJahoda 

Dunaliella acidophila 
(Kalina) Massyuk 

Subgenus Dunaliella 

Section Tertiolecta (Marine species. Optimum salinity < 6% NaCI) 

Dunaliella tertiolecta 
Butclier 



Nepal 
Nepal 

CCAP 19/35 Freshwater; acidic 

sulphurous pool, Pisciarelli, 
Naples, Italy 



CCAP19/6B Brackish; Oslo Fjord, Norway 
CCAP19/6B Brackish; Oslo Fjord, Norway 



Dunaliella quartolecta 
Butcher 

Dunaliella primolecta 
Butcher 



Dunaliella tertiolecta [this study] 
Dunaliella tertiolecta [this study] 



Dunaliella tertiolecta [this study] 

Dunaliella polymorpha Dunaliella tertiolecta [this study] 
Butcher 



Dunaliella tertiolecta [this study] 



Dunaliella maritima 
Butcher 

Section Dunaliella (Halophilic species. Optimum salinity > 

Dunaliella salina 
Teodoresco 



CCAP 19/27 
CCAP 19/27 
Dtsi 

UTEX999 

CCMP1320 

CCMP1302 

CCMP364 

FHL 

DCCBC5 
SAG 13.86 
ATCC30929 
DCCBC26 
CCAP 19/8 

UTEXIOOO 

UTEXIOOO 

hd8 

CCAP19/7C 
SAG42.89 



Italy: Venezia 
Norway: Oslofjord 
Salt flat. USA? 
Salt flat. USA? 
Salt flat. USA? 



Norway: Oslofjord 

United Kingdom: Plymouth 

Marine; Southampton, 
Hampshire, England 

English Channel, Plymouth, 
Devon, England 

English Channel, Plymouth, 
Devon, England 

China? 

Brackish; River Crouch, Essex, 
England 



I NaCI. Accumulates carotenes) 

CCAP19/18 Hypersaline; Hypersaline 
brines, Hutt Lagoon, 
Western Australia 

CCAP19/18 Hypersaline; Hypersaline 
brines, Hutt Lagoon, 
Western Australia 



CCAP 19/25 
UTEX1644 

CONC006 

CONCOOl 

CONC007 

DCCBCl 

DCCBC2 

hd6 



Point Colorado Salinas; La 
Paz, Baja California, Mexico 

Salar de Atacama, Chile 

Laguna La Rinconada, Chile 

Salar de Atacama, Chile 

Lake Tyrell, Victoria, Australia 

South Korea 

Israel 

Yucatan, Mexico 



AF3 13445 



Albertano 
(1981) 



Foyn (1928 
or earlier) 

Foyn (1928 
or earlier) 



Butcher 
(1953) 

Gross (1936) 



AY572957 

EF473748 

AY654300 

EF473730 

AF3 13435 

AF3 13433 

DQ377096 

DQ377097 

DSU66956 

AY686684 

EF473738 

EF473742 

DQ224338 

DQl 57054 

AY582942 



Gross (1936) DQ377092 



Butcher 
(1954) 



Kaethner 
(1982) 

Kaethner 
(1982) 



Loeblich 

(1967) 

(1990) 

(1990) 
Polle 



DQl 16745 
DQl 57053 

AY582086 



AF546098 

EF473746 

HMl 40783 
AF3 13429 

AF3 13425 
AF546092 
AF3 13427 
AY549442 
AY5 12973 
DQl 16743 
AF546094 



16596847 



DQ377089 87047580 
HM060646 



HM243579* 



47933783 

145587830 

55979209 

145587823 

16596842 

16596841 

87047587 

87047588 

2627284 

56578596 

145587825 

145587827 

77955899 

77539932 

50952902 

87047583 

71482600 
77539931 

51035303 

33333776 
145587829 



16596839 

16596837 
33333770 
16596838 
47499297 
46250926 
71482598 
33333772 
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Table 1 Taxonomic classification. Culture Collection, Geographic Origin and GenBank acession numbers of the strains 
included in this study. (Continued) 



Dunoliello viridis [this study] 



Dunoliello tertiolecto [this study] 

Dunaliella tertiolecto [this study] 

Dunoliello tertiolecto [this study] 

Dunoliello tertiolecto [this study] 

Dunoliello viridis [ 9, 1 0, 27] 

Dunoliello viridis [ 9, 1 0, 27] 



AC 144 

184.80 

OUC66 
"hd4" 

OUC38 
"hd3" 

OUC36 
"hd2" 

0UC21 
"hdl" 

9802 

"hd5" 

DS18S1 

DS18S2 

DS18S3 

Dsge 

CCAP 19/3 

UTEX200 

MSI-1 
ITC5100 



ITC5101 



Tanggu, China 
Tunisia, North Africa 

China 
China 
China 
China 
China ? 

Inner-Mongolia 

Mexico? 

Mexico? 

Mexico? 

Belgium: Gent 

Brackish; dirty salt lake, 
Soviet Union 

Brackish; dirty salt lake, 
Soviet Union 

Vargas, Gran Canaria, Spain 
Punta, Gran Canaria, Spain 



ITC5102 Tenefe, Gran Canaria, Spain 



ITC5103 Rio, Lanzarote, Spain 



ITC5104 Guatiza, Lanzarote, Spain 



ITC5105 Janubio, Lanzarote, Spain 



ITC5106 Carmen (Majo), 

Fuerteventura, Spain 

ITC5107 Ahana, Alava, Spain 



ITC5118 lie de Re (01), France 



ITC5122 lie de Re (05), France 



ITC51 14 La Tapa, Cadiz, Spain 



ITC5119 lie de Re (02), France 



(2005) ? 

(2005) ? 

(2005) ? 

(2005) ? 

(2007) ? 
(2005) 



AF546096 
AY549441 
AY577766 
DQl 16741 



Mainx 
Mainx 



EF695405 

DQl 16742 

FJ360756 

FJ360757 

FJ360758 

EF473732 

EF473744 



de la Jara ( 
Mendoza 
(2005) 

de la Jara & HM035354* 

Mendoza 

(2005) 

de la Jara & HM035355* 

Mendoza 

(2005) 

de la Jara & HM035356* 

Mendoza 

(2005) 

de la Jara & HM035357* 

Mendoza 

(2005) 

de la Jara & HM035346* 

Mendoza 

(2005) 

Mendoza & HM035358^ 

Trujillano 

(2003) 

de la Jara & HM035359^ 

Mendoza 

(2005) 

Carmona & HM035348* 

Mendoza 

(2006) 

Carmona & HM035347^ 

Mendoza 

(2006) 

de la Jara & HM035350* 

Mendoza 

(2007) 

Carmona & HM035349* 

Mendoza 

(2006) 



33333774 
47499296 
49073091 
71482596 



DQl 16740 71482595 



DQl 16739 71482594 



DQl 16738 71482593 



151573027 

71482597 

213958821 

213958822 

213958823 

145587824 

145587828 



AF3 13423 16596836 

GQ337903 254838316 
HM035353* 
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Table 1 Taxonomic classification. Culture Collection, Geographic Origin and GenBank acession numbers of the strains 
included in this study. (Continued) 



Aliquot of Dunaliella salina BCA421 ITC5003 Tenefe, Gran Canaria, Spain 



Dunaliella bardawil 
nomen nudum Ben- 
Amotz & Avron 



Dunaliella salina CCAP 19/30 [27], 
obtained from Dr. Joao Varela (Faro, 
Portugal) 

Dunaliella salina [27]. Reinstated 
from SAG on April 1996 

Dunaliella salina [27] 



Dunaliella salina [27] 



Dunaliella parva Lerche 



Dunaliella salina [27], Dunaliella 
tertiolecta [this study] 

Dunaliella tertiolecta [this study] 

Dunaliella viridis [9-11] 

Dunaliella tertiolecta []0,]]]Dunaliella 
quartolecta [27] 

Dunaliella tertiolecta [10,11] 

Dunaliella maritima [27], Dunaliella 
viridis [this study] 

Dunaliella tertiolecta [this study] 



ITC5000 Marine; salt pond, near 

Bardawil lagoon. North Sinai, 
Israel 

CCAP 19/30 Marine; salt pond near 

Bardawil lagoon. North Sinai, 
Israel 

ATCC30861 Marine; salt pond near 

Bardawil lagoon. North Sinai, 
Israel 

UTEX2538 Marine; salt pond near 

Bardawil lagoon. North Sinai, 
Israel 

SAG42.88 Marine; salt pond, near. 

Bardawil lagoon. North Sinai, 
Israel 

hd7 China? 

UTEX1983 Dead Sea 

CCAP 19/9 Brackish; salt marsh, Northey 
Island, Essex, England 

CCMP362 

SAG19-1 Marine: Lacul Sarat, Romania 



Mendoza 
(1992) 

Ben-Amotz 
& Avron 
(1976). 

Ben-Amotz 
& Avron 
(1976). 

Ben-Amotz 
& Avron 
(1976) 

Ben-Amotz 
& Avron 
(1978) 

Ben-Amotz 
& Avron 
(1976) 

(1973) 

Butcher 
(1956) 

Gold 

Lerche 
(Before 
1938) 



hd9 China? 
CONCOl 0 Salar de Atacama, Chile 



HM035352* 
HM03535r 

EU932917 

AF3 13431 

DQ377085 

EF473741 

DQ 11 6744 
AF313441 
AF3 13439 

AF313437 
DQ377091 

DQl 16746 
AF3 13421 



205361369 



16596840 



87047576 



71482599 
16596845 
16596844 

16596843 
87047582 



Dunaliella pseudosalina 
Massyuk & Radchenko 

Section Viridis (Halophilic species. Optimum salinity > 6% NaCI. Cells always green. Do not accumulate carotenes. Cells radially symmetrical) 



71482601 
16596835 



Dunaliella minuta Lerche Dunaliella tertiolecta [this study] 



Dunaliella bioculata 
Butcher 



Dunaliella viridis Massyuk 



Dunaliella tertiolecta [this study] 
Dunaliella tertiolecta [this study] 

Dunaliella tertiolecta [this study] 



CCAP 19/5 

SAG23.86 
CCAP 19/4 

UTEX199 

UTEX199 

CONC002 
CONC002 
SAG44.89 
SHU 

DCCBC4 
DCCBC3 



Marine; sand and sea water, Jowett 
Roscoff, France (1967) 

Brackish; salt lake, Soviet Mainx 
Union 

Brackish; salt lake, Soviet Mainx 
Union 

Brackish; salt lake, Soviet Mainx 
Union 

Salar de Atacama, Chile (1990) 
Salar de Atacama, Chile (1990) 

China? 

Great Salt Lake, Utah, USA 
Great Salt Lake, Utah, USA 



HM035345* 

AY582085 
HM035344* 

DQl 57433 

DQ377086 

AF313419 
DQ377098 
87047600 
AY878700 
AY686685 
AY828227 



51035302 



76097092 

87047577 

16596834 
87047589 
61200914 
58339343 
56578597 
61200913 



Section Peirceinae ((Halophilic species. Optimum salinity > i 

Dunaliella percei Nicolai Dunaliella tertiolecta [10,1 1] 
& Baas-Becking 

Dunaliella tertiolecta [10,11] 

Unknown Dunaliella Species 
Dunaliella sp. Dunaliella tertiolecta [27] 

Dunaliella salina [27] 



Dunaliella tertiolecta [this study] 
Dunaliella tertiolecta [this study] 



NaCI. Cells always green. Do not accumulate carotenes. Cells bilaterally symmetrical) 



CCAP 19/2 Brackish; California, USA 



UTEX21 92 Brackish; California, USA 



CCAP19/23 Marine; 

CCAP19/12 Brackish; North Sinai, Israel 

CCMP367 Salt flat 

CCMP220 Salt flat 

CCMP1923 Salt flat 

CCMP1641 Salt flat 



Nicolai 
(1931) 

Nicolai 
(1931) 

Pennick 

Ginzburg 
(1976) 



HM035343* 



AF3 13443 



HM03534r 
HM035342^ 

DQ377087 
DQ377095 
DQ377094 
DQ377093 



16596846 



87047578 
87047586 
87047585 
87047584 
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Table 1 Taxonomic classification. Culture Collection, Geographic Origin and GenBank acession numbers of the strains 
included in this study. (Continued) 



Dunoliello tertiolecta [this study] 

Dunaliella viridis [this study] 
Dunaliella viridis [this study] 
Dunaliella viridis [this study] 
Dunaliella salina [this study] 

Dunaliella tertiolecta [this study] 
Dunaliella viridis [this study] 

Dunaliella salina [this study] 
Dunaliella viridis [this study] 



SAG 19.6 

FLl 

BSFl 

BSF2 

BSF3 

006 



hdlO 

ABRIINW 
Ml/1 

ABRIINW 
Ml/2 

ABRIINW 
Ul/1 

ABRIINW 
U2/1 

SPM0112- 
3 

SPMO201- 
3 

SPM0128- 
2 

SPMO109- 
1 

SPM0112- 
4 

SPMO207- 
3 

SPMO200- 
3 

SPMO201- 
4 

SPMO201- 
5 

SPMO201- 
6 

SPM0112- 
1 

SPM0112- 
2 

SPMO300- 
4 



Salt flat 

USA: Utah, Bonneville Salt 
Flats 

USA: Utah, Bonneville Salt 
Flats 

USA: Utah, Bonneville Salt 
Flats 



William 
Henley 

William 
Henley 

William 
Henley 

A.W. 

Coleman, U. 
Brown 



Iran? 
Iran? 
Iran? 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 

Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 



AY582086 
DQ377099 
DQ377081 



AF033278 

DQ 116747 
EU927374 

EU927373 

FJ 164063 

FJ 164064 

DQ377103 



DQ377115 



DQ377111 



DQ377112 



DQ377113 



51035303 
87047590 
87047572 



DQ377082 87047573 



DQ377083 87047574 



2645739 

71482602 
197290927 

197290646 

205371718 

205371719 

87047594 



DQ377110 87047601 



DQ377105 87047596 



DQ377105 87047596 



DQ377104 87047595 



87047606 



DQ377107 87047598 



87047602 



87047603 



87047604 



DQ377101 87047592 



DQ377102 87047593 



DQ377118 87047609 
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Table 1 Taxonomic classification. Culture Collection, Geographic Origin and GenBank acession numbers of the strains 
included in this study. (Continued) 







SPM02in- 

wJI 1 V 1 vyz. 1 \J 

3 


Salt flat USA- Oklahoma Salt 
Plains National Wildlife 
Refuge 




□03771 16 


87047607 

kJ I / \J\J / 




DunoliellQ vifidis [this study] 


SPMO200- 
8 


Salt flat USA- Oklahoma Salt 
Plains National Wildlife 
Refuge 




□0377108 


87047S99 




Diinnliplln viririi'^ fthi'^ '^tiiHvl 

Ly L///Of//C//0( V 1 1 IKJ IS) l_LI 1 IS) J LU KJy J 


1 


Salt flat USA- Oklahoma Salt 
Plains National Wildlife 
Refuge 




□03771 21 

L-/ / / 1 Z_ 1 


«7n4761 2 




Dunoliello vifidis [this study] 


SPMO200- 
2 


Salt flat USA- Oklahoma Salt 
Plains National Wildlife 
Refuge 




□03771 06 


87047597 




Dunciliellci vifidis [this study] 


SPMO201- 
2 


Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 




DQ377109 


87047600 




Dunoliello vifidis [this study] 


SPMO202- 
4 


Salt flat USA- Oklahoma Salt 
Plains National Wildlife 
Refuge 




DQ3771 14 


87047605 




Dunoliello vifidis [this study] 


SPMO300- 
5 


Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 




DQ3771 19 


87047610 




Dunoliello vifidis [this study] 


SPMO600- 
1 


Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 




DQ377120 


8704761 1 




Dunoliello vifidis [this study] 


SPMO BPS 


Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 




DQ377122 


87047613 




Dunoliello vifidis [this study] 


SPMO 
980625-IE 


Salt flat, USA: Oklahoma, Salt 
Plains National Wildlife 
Refuge 




DQ377123 


87047614 


vJl\ Icl yi uupb 














L./ ilUi I lyUUi 1 lUi IUj 

reinhardtii 




SAGl 1-31 








1 1 1 07^40S 


Chlamydomonas 
allensworthii 




isolate 7 






AF326855 


13274452 


Oedogonium nodulosum 










DQ078301 


71482662 


Bulbochoete rectonguloris 
var. hiloensis 




UTEX954 


Catawmont, Massachusetts, 
USA 


Cook 1962 


AY962677 


62183531 


Astrephomene 
gubernoculifero 




UTEX2479 






AGU66932 


2982755 


Pandorina morum 




EP3 






AF378359 


14165185 


Eudorina elegons 




UTEX 1192 






AF098173 


4007503 


Yamagishiella unicocca 




UTEX2431 






AF375785 


14091689 


Phocotus lenticuloris 




970805-20 






AY009933 


18025134 


Lobochlomys cu Ileus 




UTEX1057 


Maxville, Florida, USA 


Smith 


CCU66946 


2627275 



Asterisks indicate the strains sequenced in tliis study. Underlined strains correspond to identical isolates stored in different culture collections. 



Acronyms: Culture Collection of Algae and Protozoa, UK (CCAP); Sammlug von Algenkulturen, Germany (SAG); University of Texas Culture Collection of Algae, USA 
(UTEX); American Type Culture Collection, USA (ATCC); Dunaliella Culture Collection at Brooklyn College, USA (DCCBC); Provasoli-Guillard National Centre for the 
culture of Marine Phytoplankton, USA (CCMP); Universidad de Concepcion, Chile (CONC); Banco Canario de Algas, Spain (BCA) ; Instituto Tecnologico de Canarias, 
Spain (ITC). 



not able to elucidate completely the species boundaries 
within the different groups of Dunaliella, since in some 
cases it was observed that there was a lack of CBCs 
between known distinct species. 

In the special case of D. salina, high variation levels 
have been reported for decades [19-24]. However, only 
Massjuk [18] translated geographical, physiological, and 



morphological variables into the recognition of two sub- 
species (D. salina sp. salina and D, salina sp. sibirica 
Massjuk and Radch.) and three forms {D, salina sp. sal- 
ina f. salina, D. salina sp. salina f. oblonga Lerche, and 
D. salina sp. salina f. magna Lerche). Later on, ITS2 
phylogenetic analyses suggested the existence of two dis- 
tinct phylogenetic species within the taxonomic entity 
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currently known as D. salina [10,11,32], indicating the 
possibility of cryptic speciation [10]. Our ITS2 phyloge- 
netic analysis does confirm the existence of three differ- 
ent groups within D. salina; however, the CBC results 
did not resolve if these groups may correspond to dis- 
tinct species, although several strains of each D. salina 
group shared one CBC with D, salina strains in other 
groups. On the other hand, the high morphological and 
physiological variability found within the Spanish D. sal- 
ina strains under the same lab conditions [24] was not 
correlated with the phylogenetic observations of this 
study. This finding indicates that the phylogenetic close- 
ness found with ITS sequences does not reflect common 
physiological or morphological attributes. Moreover, our 
data unequivocally suggest that D. salina is not mono- 
phyletic, at odds with previous hypotheses [10,11,32]. 

The other objective of this work was to elucidate if the 
physiological uniqueness found in "Janubio" [such as its 
unique fatty acid profile and accumulation of high levels 
of carotenes under low light flux density conditions; 
Mendoza et al.: A new strategy for carotenogenesis 
under conditions of cellular stress in Dunaliella (a 
Potential New Species), submitted] could be confirmed 
by the ITS2 data and the CBC analysis. However, 
although we observed that this strain has a unique ITS2 
sequence profile, and had more than one CBC with the 
other phylogenetically related D, salina^ our data do not 
allow us to conclude that this is a new species, and 
further studies must be performed to find out if the dif- 
ferences observed are just reflecting a high intra- specific 
variability. Finally, in agreement with previous studies 
[8,11,25], our ITS2 data failed to furnish evidence for 
isolation by distance among D, salina strains. 

Conclusion 

This work demonstrates that the taxonomy of Duna- 
liella should be revised. The great diversity observed 
within the ITS2 sequences of D. salina suggests that dif- 
ferent biological groups may exist within this taxon; 
however, this was not confirmed with the CBC analysis. 
Likewise, although the Spanish D. salina strain ITC5105 
"Janubio" was characterized by a unique ITS2 sequence, 
the hypothesis that it may be a new species could not 
be confirmed by the CBCs analysis, requiring further 
morpho-physiological and genetic investigation. Overall, 
the use of CBCs to define species boundaries within 
Dunaliella was not conclusive in some of the cases 
assessed. 

Methods 

Strains, DNA extraction and ITS amplification 

We sequenced the ITS2 region of 13 D. salina strains 
from Spanish and French saltworks, one strain obtained 
from the Culture Collection of Algae and Protozoa UK 



(CCAP), and one D. salina strain that has been main- 
tained in the Instituto Tecnologico de Canarias (ITC 
henceforth) for several years (purchased from CCAP as 
Dunaliella salina 19/30). We also sequenced other 
Dunaliella species {D. minuta CCAP 19/5, D. tertiolecta 
CCAP19/23 and CCAP19/6B, D. bioculata CCAP19/4) 
[Table 1]. The sequences of the other strains analyzed 
in this study were obtained from the ITS2 Database 
(http://its2.bioapps.biozentrum.uni-wuerzburg.de/? 
about). Detailed information about the strains used in 
this study can be found in Table 1. 

DNA extraction was performed with a chelex-100 
(Biorad, CA, USA) resin-based protocol [33]. For the 
DNA amplification of the ITS region, primers AB28 and 
TW81 in Goff et al. (1994) [34] were used. DNA ampli- 
fication was carried out in a total volume of 25 [il with 
IX iQ SYBR Green Supermix (Biorad, CA, USA) and 10 
pM of each primer in a Smart Cycler thermocyler (Cep- 
heid,CA, USA) as follows: 5 min at 94°C; 5 cycles of 1 
min at 94°C, 2 min at 50°C and 1 min at 72°C; 30 cycles 
of 1 min at 94°C, 1 min at 62°C and 1 min at 72°C, with 
a final extension of 5 min at 72°C. 

PCR products were first electrophoresed in a 1.5% 
agarose gel to assure that a single band of 500-600 bp 
was present, then purified using the Real Clean Spin kit 
(REAL, Durviz S.L.U., Valencia, Spain), and finally bi- 
directionally sequenced on an ABI PRISM 3730x1 auto- 
matic sequencer (Applied Biosystems, CA, USA) at the 
DNA sequencing services of Macrogen (Korea). 

Phylogenetic analyses 

Sequences and their individual secondary structures 
were obtained from the ITS2 Database [35-37]. Newly 
obtained ITS2 sequences were annotated according to 
Keller et al. [38], and their secondary structures pre- 
dicted by homology modeling [39]. The phylogenetic 
analysis followed the procedure outlined in Schultz and 
Wolf [7] in accordance with Keller et al. [5]. The soft- 
ware used for the ITS2 sequence-structure analysis can 
be obtained from http://its2.bioapps.biozentrum.uni- 
wuerzburg.de/?about. A global, multiple sequence-struc- 
ture alignment was generated in 4SALE vl.5 [40,41]. 
The sequences and their individual secondary structures 
were synchronously aligned making use of an ITS2 
sequence-structure specific scoring matrix [40], and the 
start and end of the alignment was manually adjusted. 
Based on primary and secondary structure information, 
phylogenetic relationships were reconstructed by Prof- 
DistS, through the use of an ITS2 specific, general time 
reversible substitution model [42,43]. Bootstrap support 
[44] was estimated on 100 pseudo-replicates. The result- 
ing tree was visualized with TreeView [45]. 

To study the species boundaries within Dunaliella we 
followed the "distinguishing species" instructions [6] 
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based on compensatory base changes (CBCs) in the 
ITS2 secondary structure, and we used the CBCAnalyzer 
option implemented in 4SALE. 

Additional material 



Additional file 1: Compensatory Base Changes (CBC) analysis. 

Compensatory Base Changes (CBCs) between different groups and 
species witlnin the Dunaliella taxon (Excel file). 
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